Reorganization energy of electron transfer at the solvent glass transition 
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We present a molecular-dynamics study of the solvent reorganization energy of electron transfer 
in supercooled water. We observe a sharp decrease of the reorganization energy at a temperature 
identified as the temperature of structural arrest due to cage effect as discussed by the mode cou- 
pling theory. Both the heat capacity and dielectric susceptibility of the pure water show sharp 
drops at about the same temperature. This temperature also marks the onset of the enhancement 
of translational diffusion relative to rotational relaxation signaling the break-down of the Stokes- 
Einstein relation. The change in the reorganization energy at the transition temperature reflects 
the dynamical arrest of the slow, collective relaxation of the solvent related to Debye relaxation of 
the solvent dipolar polarization. 



I. INTRODUCTION 

Theories of activated chemical dynamics and trans- 
port phenomena in condensed phase are often based 
on transition-state ideas invoking equilibrium thermody- 
namics to describe the reaction flux across the transition- 
state surface separating the reactants from the products. 
The Marcus-Hush theory of electron transfer (ET) reac- 
tions fully relies on the transition-state formalism defin- 
ing the ET activation barrier in terms of two thermody- 
namic parameters, the free energy gap and the nuclear 
reorganization energy. The former is the difference in 
free energies between the final and initial ET states, and 
the latter determines the curvature of two free energy 
parabolasi When the donor-acceptor energy gap AE is 
chosen as the reaction coordinate, the classical reorga- 
nization energy is given through the energy gap second 
cumulant: 



A, = {{6AEf)/2kBT. 



(1) 



The transition-state description becomes inapplica- 
ble when the time of passage of the activation barrier 
is comparable to the relaxation time of the condensed 
medium (solvent). The population of the activated state 
gets depleted and one arrives at the friction-affected 
chemical kinetics described by the Kramers theorji^ and 
its modifications!^ For ET, this regime corresponds to 
solvent-controlled reactions when the preexponential fac- 
tor of the rate is inversely proportional to a solvent relax- 
ation time.'* One can anticipate the next step in this hi- 
erarchy of relaxation times when the time of the reaction 
itself (not just the time of barrier passage) becomes com- 
parable to the solvent relaxation time. Such conditions, 
which apply to ultra-fast reactions in high-temperature 
solvents and to reactions in slowly-relaxing viscous sol- 
vents, will result in the loss of ergodicity of the system 
and the break-down of the equilibrium description of the 
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reaction activation barrier (in contrast to the alteration 
of the rate preexponent in the Kramers description).^ A 
general theory of chemical rates at such conditions is still 
missing even though nonergodic behavior may apply to 
a broad class of reactions in supercooled liquids and in 
biopolymersii& The latter case is particularly relevant to 
the problem of nonergodic activation since the dynamics 
of biopolymers is characterized by a broad spectrum of 
relaxation times, ^'^ and, for a given reaction rate, at least 
a subset of nuclear modes may become nonergodic. 

As the second cumulant of the difference in the solute- 
solvent interaction potential (eq^, the solvent reorgani- 
zation energy of ET is related to a response function of 
the solvent to the presence of the solute. The latter can 
be related to the susceptibility of the pure solvent cor- 
responding to a nuclear collective coordinate coupled to 
the solute electronic states (dipolar polarization and den- 
sity for dipolar solvents, quadrupolar polarization and 
density for non-dipolar solvents) Susceptibilities, which 
are second derivatives of a thermodynamic potential, are 
known to show specific behavior at points of thermody- 
namic instability (phase transitions) or at the onset of 
nonergodicity (glass transition). It is of general inter- 
est to understand how chemical reactions are affected by 
these special points. For glass transitions, the heat ca- 
pacity of a glassformer (the second cumulant of energy 
at constant volume or the second cumulant of entropy at 
constant pressure) passes through a peak which sharp- 
ens with increasing fragility of the liquid^ In this Letter 
we present simulation results indicating that the solvent 
reorganization energy of ET demonstrates a similar be- 
havior. It first increases with lowering temperature and 
then sharply drops at the temperature of kinetic glass 
transition. 



II. SIMULATION RESULTS 

We perform NVE and NPT molecular-dynamics 
(MD) simulations for a system composed of one solute (p- 
nitroaniline) and N — 466 water molecules, with periodic 
boundary conditions. The extended simple point charge 
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model (SPC/E) is adopted for water, p-nitroaniline is 
a small charge-transfer dye with a charge-transfer elec- 
tronic transition resulting in about 3.7 D change in the 
dipole moment^ Its interaction with water is a sum of 
Coulomb interactions between solute and solvent par- 
tial charges and Lennard-Jones (LJ) potentials for which 
Lorentz-Berthelot combination rules are used. The NVE 
simulations are done for a cubic simulation cell with 
the side length of 24.075 A and water density of 0.997 
gm/cm^. The cutoff radius of 12.0 A was employed. 
Simulations have been carried out in the range of tem- 
peratures from supercooled region (50 K) to superheated 
region (509 K) liSiiiiii^ In parallel to solvation simula- 
tions, simulations of the pure SPC/E water {N = 466) 
at the same thermodynamic conditions have been done. 
Equilibration has been carried out over a time period of 
400 ps and production runs have been accumulated over 
the duration of 1.0 ns at 509 K up to 6.0 ns at 50 K from 
configurations stored at the interval of 0.1 ps. 
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FIG. 1: Reorganization energy (a), heat capacity (b), and 
dielectric susceptibility (c) vs T. In (a), circles refer to NVE 
simulations at p = 0.9997 cm'^/g, squares indicate the NPT 
simulations at P = 1 atm, diamonds refer to the fast Gaussian 
component of the solvent reorganization energy, Xq- Dotted 
lines are drawn to guide the eye. 

Atomic charges of p-nitroaniline were obtained by 
fitting the electrostatic potential from ab initio elec- 
tron structure calculations using GAUSSIAN'03i& The 
ground state geometry of p-nitroaniline was obtained on 
the MP2 level (6-3H-G*) using X-ray data^^ for the ini- 
tial geometry. The excited-state charge distribution is 
obtained from SCI calculations. The ground, 7.18 D, 
and excited, 10.88 D, dipolcs are in reasonable agreement 
with the literature dataiii*i^ Coordinates and charges of 
p-nitroaniline in the ground and excited states are listed 
in the supplement. 

Figure Q^a) shows the solvent reorganization energy 



calculated according to eq ^ in which the energy gap 
AE is replaced with the difference of the solute-solvent 
interaction potential in the charge-transfer (5*4) and 
ground {So) electronic statesiii The average is performed 
over the simulation trajectories generated for the so- 
lute in the ground state. Circles indicate constant- 
volume simulations and squares refer to simulations at 
constant pressure of P = 1 atm. The reorganization 
energy increases with lowering temperature in the high- 
temperature liquid^ and then turns down and drops 
sharply at T* ~ 219 K. 
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FIG. 2: Stokes shift correlation function (eqllj of ground-state 
p-nitroaniline in SPC/E water. 

Shown in Figure mb,c) are also the heat capacity, 

cv/kB = 3/2 + {{5Ef)lN{kBTf (2) 

and dielectric susceptibility 

XP^{kBTV)-'{{5mf) (3) 

of SPC/E water {E is the total energy and M is the total 
dipole moment of a sample of liquid with volume V con- 
taining N molecules). Both solvent susceptibilities show 
the dependence on temperature very similar to that of As . 
Figure n thus stresses the similarity in the temperature 
variation of the ET reorganization energy and suscepti- 
bilities of the pure solvent at the onset of nonergodicity 
at T = T*. 

The reorganization energy does not drop to zero, 
but instead levels off at low temperatures. The low- 
temperature component of As is related to fast phonon- 
like solvent modes which do not become dynamically ar- 
rested below T* . The fast reorganization component is 
clearly seen as the initial Gaussian decay of the Stokes 
shift correlation function^^ (Figure 

S{t)^C{t)IC{Q), C{t) = {AE{t)AE{Q)). (4) 

S{t) was fitted to a biphasic form containing the Gaussian 
(G) and stretching exponential (E) parts 

S{t) = ^Ge"(*/"°)' + (1 - ^G)e~(*/"^^' . (5) 

The stretching exponent (5 obtained from the fit is equal 
to one above 272 K and starts to drop below this tem- 
perature reaching the value of 0.34 at 219 K (see the 
supplement). 
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From the fit of the shiiulated Stokes shift data, the 
Gaussian relaxation timeiS tq — 70 fs is weakly temper- 
ature dependent and is approximately given by a linear 
function of T (K): ln(TG(T)/ps) = -1.34 - 0.0042 x T 
(Figure |2Jl. In contrast, the exponential relaxation time 
te increases sharply with lowering T and can be approx- 
imated by the Vogel-Fulcher (VF) law: ln(T£;(T)/ps) = 
-3.73 + 447/(r - 161), where T is in K3 The Gaus- 
sian reorganization component, Ag — AqXs, turns out 
to be weakly dependent on temperature (Figure ^a), di- 
amonds). The combination of the weak temperature de- 
pendence of Xg with the ultra- fast relaxation of the Gaus- 
sian component of S{t) clearly indicates that it is the slow 
exponential component of S(t) that becomes dynamically 
arrested at the transition to nonergodicity. The drop of 
As at r < T* is thus equal to A^ = (1 - Ag)\s- 
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FIG. 3: Upper panel: Relaxation times tq, te and ri vs T. 
The inset shows natural logarithm of the diffusion coefficient 
D (m^/s). The solid lines indicate the VF fits. Lower panel: 
The product of the self-diffusion coefficient and the rotational 
relaxation time n vs T. The inset shows the maximum value 
ce2{t*) of the non-Gaussian parameter a2{t) defined by eq|H| 

The slow component of the Stokes shift dynamics can 
be identified with collective orientational relaxation of 
the dipolar polarization of water related to the frequency- 
dependent dielectric constant e(w). The variation of e(w) 
with lowering temperature is characteristic of many polar 
glassformers (Figure^J).^^ In particular, the peak ujmax of 
the dielectric loss e"(ijj) shifts to lower frequencies with 
cooling the solvent. The width of the dielectric spectrum, 
however, does not change, and the dielectric loss data at 
different temperatures can be superimposed on one mas- 
ter curve by proper re-scaling (Figure 0] upper panel). 
The increase in the width of dielectric loss is commonly 
associated with liquid heterogeneity/^ Figure 0] thus in- 
dicates that dielectric response is homogeneous. 

The collective relaxation of solvation shells reflected by 
Te (T) can be compared to the single-particle correlation 
function of the dipole vector m(t): 

Ci(t) = (m(i)-m(0))/(m(0)2). (6) 



The relaxation time ti{T) for Ci{t) shows a super- 
Arrhenius temperature dependence similar to te{T) 
(Figure |3J which can be approximated by the VF law: 
ln(ri(T)/ps) = -2.15 + 578/(7-137). 

What is the nature of the dynamical arrest occurring 
below ~ 219 K? Previous extensive simulations of SPC/E 
water have shown a power-law temperature dependence 
of transport coefRcients with the critical temperature of 
ca. - 186 - 200 K.22.23,24,25 gg^gg^j |-j^jg observa- 
tion and scaling laws for the intermediate scattering func- 
tion, the critical temperature was assigned to the ideal 
ergodic-to-nonergodic kinetic glass transition predicted 
by the mode-coupling theory (MCT) 1^ Our data for the 
self-diffusion coefficient of SPC/E water D{T) (obtained 
from MD trajectories using the Einstein relation) em- 
brace a broader range of temperatures than those re- 
ported in refs and The self-diffusion coefficient 
(Figure E| inset) can be fitted by both the power-law 
(ln(D(T) X s/m2) = -19 -I- 2.51n(T/181 - 1)) and VF 
(ln(L»(T) X s/m2) ^ -16.5 + 540/(T - 136)) functions, 
with the latter providing a better global fit. The limiting 
VF temperatures for single-particle correlations reflected 
by D{T) and Ci(T) (136 K and 137 K, respectively) are 
significantly lower than the corresponding temperature 
(161 K) from the many-particle correlation time T£;(T)iSi 

The temperature T* at which As(T) and solvent sus- 
ceptibilities start to drop also marks the onset of the sep- 
aration of translational and rotational diffusion signaling 
the breakdown of the Stokes-Einstein relation (Figure 01 
lower panel) 

i:'(T) X Ti(T) ~ const (7) 

Indeed, the product of the diffusion coefficient and the ro- 
tational relaxation time is approximately constant down 
to the temperature ~ 207 K.^^ If 165 K is adopted 
as the glass transition temperature for water^ the on- 
set of translational enhancement falls in the range 1.2- 
1.3 Tg found in laboratory experiment. This effect is 
commonly explained by either static spatial^''''^"'^^''^^ or 
dynamioS heterogeneity. In the frustration-limited do- 
main picture of Tarjus and Kivelson^SSiffi the turning 
temperatures in Figures ^ a-^d El can be associated with 
the onset of domain formation. This interpretation is 
questionable, however, given the small size of our sim- 
ulation box which cannot incorporate mesoscopic do- 
mains. Note also that no discontinuous change in pair 
distribution functions is observed at T* except previously 
observedSSi^ sharpening of solvation shell peaks. 

The results of our simulations better fit the picture 
of spatially heterogeneous dynamics which assumes the 
presence in a supercooled liquid of groups of mobile 
molecules 1^2*2^4^ These groups may represent clusters of 
hydrogen bond defects in SPC/E water, chains of mo- 
bile particles in monoatomic LJ fiuids}^ or some other 
structures. The collective motion of such mobile hetero- 
geneities provides the enhancement of translational dif- 
fusion which is correlated with the increase of the non- 
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FIG. 4: Real (e'(tj)) and imaginary (e"(a;)) parts of the 
frequency-dependent dielectric constant of SPC/E water at 
different temperatures. The upper panel shows the re-scaled 
dielectric loss e"((jj)/e(a;max) vs logiQ(tj/aJmax) calculated at 
198, 219, 242, and 300 K; ojniax is the frequency maximum 
of the dielectric loss function. The data points at different 
temperatures cannot be distinguished on the plot scale. 



lations. On the other hand, dielectric loss data show no 
heterogeneity suggesting a more homogeneous distribu- 
tion of molecular rotations. 

This interpretation is consistent with recent dielectric 
data by Richert et alm^ which do not show any increase 
in spatial heterogeneity in the region of temperatures 
Tg < T < l-2Tg where previous reports^ had indicated 
the breakdown of the Stokes-Einstein relation. Although 
no broadly accepted explanation of this phenomenon cur- 
rently exists, the onset of translation/rotation decoupling 
is normally associated with the critical MOT tempera- 
ture. We will therefore resort to the interpretation of 
temperature T* at which As dips to its Gaussian com- 
ponent as the point of kinetic transition to nonergodicity 
with many features of the critical temperature of an ideal 
glass transition predicted by MCTi^ 



III. CONCLUDING REMARKS 

The MCT critical temperature Tc located above the 
calorimetric glass transition temperature marks the 
change in the mechanism of relaxation in viscous liquids. 
Above Tc, transport phenomena and response to a solute 
field are dominated by the rattling and librations of par- 
ticles in self-consistently maintained cages. Below Tc, on 
the other hand, the particles are almost arrested in the 
free-energy landscape and transport and relaxation are 
triggered by thermally activated hopping and rotation 
over saddle points. The hopping mechanism is recorded 
by experiments with the observation window up to 10^ 
s in the temperature range Tg < T < Tc. Many redox 
reaction occur on the nanosecond and faster time-scale. 
Since hopping relaxation does not occur on the reaction 
time-scale, MCT formalism may appear to be appropri- 
ate for describing ET reactions in viscous supercooled 
liquids. The sharp drop of the solvent reorganization 
energy close to Tc found in our simulations is therefore 
expected to mimic the laboratory ET experiment. 



Gaussian parameter— 

a2(t) = 3(r4(t))/5(r2(i))2-l, (8) 

where is the mean square displacement. The max- 

imum value a2{t*) (Figure|311ower panel inset) is reached 
at time t* . Both a2{t*) and t* grow with cooling indicat- 
ing increasing dynamic heterogeneity of molecular trans- 



Acknowledgments 

This research was supported by the National Science 
Foundation (CHE-0304694). This is pubhcation #643 
from the ASU Photosynthesis Center. Useful discussions 
with Prof. C. A. Angell and Prof. R. Richert are grate- 
fully acknowledged. 



^ Marcus, R. A. Rev. Mod. Phys. 1993, 65, 599. 
^ Gardiner, C. W. Handbook of stochastic methods; Springer: 
Berlin, 1997. 

^ Grote, R. F.; Hynes, J. T. J. Chem. Phys. 1980, 73, 2715. 
* Barzykin, A. V.; Frantsuzov, P. A.; Seki, K.; Tachiya, M. 



Adv. Chem. Phys. 2002, 123, 511. 
^ Matyushov, D. V. J. Chem. Phys. 2005, 122, 084507. 
^ Dick, L. A.; Malfant, I.; Kuila, D.; Nebolsky, S.; No- 

cek, J. M.; Hoffman, B. M.; Ratner, M. A. J. Am. Chem. 

Sac. 1998, 120, 11401. 



5 



Frauenfelder, H.; Parak, F.; Young, R. D. Ann. Rev. Bio- 

phys. Biophys. Chem. 1998, 17, 451. 
* Andreatta, D.; Perez, J. L.; Kovalenko, S. A.; Ernst- 

ing, N. P.; Murphy, C. J.; Coleman, R. S.; Berg, M. A. 

J. Am. Chem. Soc. 2005, 127, 7270. 
^ Matyushov, D. V. J. Chem. Phys. 2004, 120, 7532. 
1° Angell, C. A. Science 1995, 267, 1924. 

Farztdinov, V. M.; Schanz, R.; Kovalenko, S. A.; Ernst- 

ing, N. P. J. Phys. Chem. A 2000, 104, 11486. 

Stanley, H. E.; Buldyrev, S. V.; Canpolat, M.; 

Mishima, O.; Sadr-Lahijany, M. R.; Scala, A.; 

Starr, F. W. Phys. Chem. Chem. Phys. 2000, 2, 1551. 
" Debenedetti, P. G.; Stanley, H. E. Phys. Today 2003, 56, 

40. 

Vega, C; Sanz, E.; Abascal, J. L. F. J. Chem. Phys. 
2005, 122, 114507. 

Giovambattista, N.; Angell, C. A.; Sciortino, F.; Stan- 
ley, H. E. Phys. Rev. Lett. 2004, 93, 047801. 
Frisch, M. J.; Trucks, G. W.; Scuseria, G. E.; 
Robb, M. A.; Cheeseman, J. R.; Zakrzewski, V. G. GAUS- 
SIAN 03; Pittsburgh, PA: 2003. 

" Donohue, J. Acta Cryst. 1956, 9, 960. 

1* Rashid, A. N. J. Mol. Struc. 2004, 681, 57. 

^® Jimenez, R.; Fleming, G. R.; Kumar, P. V.; Maron- 
celh, M. Nature 1994, 369, 471. 

The Vogel-Fulcher temperature for te{T) falls slightly be- 
low Tg ~ 165 K recently identified as the glass transition 
temperature of water, see Y.-Z. Yu and C. A. Angell, Na- 
ture (London) 2004, 427, 717. 

21 Dixon, P. K.; Wu, L.; Nagel, S. R.; Williams, B. D.; 
Carini, J. P. Phys. Rev. Lett. 1990, 65, 1108. 

22 Gallo, P.; Sciortino, F.; Tartaglia, P.; Chen, S.-H. Phys. 
Rev. Lett. 1996, 76, 2730. 

2^ Sciortino, F.; Gallo, P.; Tartaglia, P.; Chen, S.-H. Phys. 



Rev. E 1996, 54, 6331. 

2* Fabbian, L.; Latz, A.; Schilling, R.; Sciortino, F.; 
Tartagha, P.; Theis, C. Phys. Rev. E 1999, 60, 5768. 

2^ Theis, C; Sciortino, F.; Latz, A.; Schilling, R.; 
Tartagha, P. Phys. Rev. E 2000, 62, 1856. 

2® Gotze, W. Aspects of structural glass transitions. In Liq- 
uids, freezing and glass transition, Vol. 1; Hansen, J. P.; 
Levesque, D.; Zinn- Justin, J., Eds.; Elsevier: Amster- 
dam, 1991. 

2^ The temperature Tc = 207 K from the plot of D{T)ri{T) 
vs T is very close to the estimate for the critical tempera- 
ture Tc ~ 208 K from the molecular MCT with a restricted 
set of correlation functions (dipole approximation) i^i This 
agreement is, however, fortuitous since the full version of 
the molecular MCT leads to a substantially higher tem- 
perature, Tc ~ 279 K.25. 

2** Chang, L; Sillescu, H. J. Phys. Chem. B 1997, 101, 8794. 

2^ Tarjus, G.; Kivelson, D. J. Chem. Phys. 1995, 103, 3071. 
Kivelson, D.; Tarjus, G. Phil Mag. B 1998, 77, 245. 

^1 Xia, X.; Wolynes, P. G. J. Phys. Chem. B 2001, 105, 
6570. 

^2 Swallen, S. F.; Bonvallet, P. A.; McMahon, R. J.; Edi- 
ger, M. D. Phys. Rev. Lett. 2003, 90, 015901. 
Jung, Y. J.; Garrahan, J. P.; Chandler, D. Phys. Rev. E 
2004, 69, 061205. 

Donati, C.; Glotzer, S. C; Poole, P. H.; Kob, W.; Plimp- 
ton, S. J. Phys. Rev. E 1999, 60, 3107. 
Schober, H. R. Phys. Chem. Chem. Phys. 2004, 6, 3654. 
Giovambattista, N.; Mazza, M. G.; Buldyrev, S. V.; 
Starr, F. W.; Stanley, H. E. J. Phys. Chem. B 2004, 
108, 6655. 

^"^ Richert, R.; Duvvuri, K.; Duong, L.-T. J. Chem. Phys. 
2003, 118, 1828. 



